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Abstract 

We analyze the chiral transport terms in relativistic superfluid hydrodynamics. In addition 
to the spontaneously broken symmetry current, we consider an arbitrary number of unbroken 



symmetries and extend the results of arXiv:1105.3733. We suggest an interpretation of some of the 
new transport coefficients in terms of chiral and gravitational anomalies. In particular, we show 
that with unbroken gauged charges in the system, one can observe a chiral electric conductivity - a 
current in a perpendicular direction to the applied electric field. We present a motivated proposal 
for the value of the associated transport coefficient, linking it to the triangle anomaly. Along the 
way we present new arguments regarding the interpretation of the anomalous transport coefficients 
in normal fluids. We propose a natural generalization of the chiral transport terms to the case of 
an arbitrary number of spontaneously broken symmetry currents. 
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I. INTRODUCTION 



The most remarkable property of liquid helium below the A-point is superfluidity. It 
is the ability of the fluid to flow inside narrow capillaries without friction, discovered by 
Kapitza [1]. The theoretical basis for understanding the phenomenon of superfluidity was 
given by Landau The hydrodynamics of a superfluid consists of two motions: the 
motion of the normal part of the fluid, and the motion of the superfluid part which is an 
irrotational one, i.e. its velocity is curl free 

as 

. A superfluid can be described fluid 
with a spontaneously broken symmetry, where the superfluid component is the condensate, 
and its velocity is proportional to the Goldstone phase gradient. The hydrodynamics of 
relativistic superfluids has been studied in |4|, and is relevant to the study of neutron stars 
5| and highly dense quark matter at the low temperature Color-Flavor locked phase 6| 
(for a general perspective on the CFL phase, see 7, sj]). A gravitational holographic dual 
description of relativistic superfluid hydrodynamics has been proposed in jl], 10]. 



Quantum anomalies in the microscopic gauge theory imply new transport terms in the 



fluid dynamics 
studied in 
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he anomalous transport in normal (i.e. not super-) fluid dynamics has been 
15||, and in the holographic gravitational framework in If]- 18]. Experimental 
signatures of this anomalous transport were proposed in [ii], Q • 

The aim of this paper is to st udy chiral effects and anomalous transport in superfluid 
hydrodynamics. In a recent work 21] , the entropic constraints on superfluid transport terms 
were analyzed, and the allowed transport terms were listed (a partial list has been obtained 



also in 



22j). This was done for the case of a single abelian spontaneously broken charge. In 



the present work, we extend the calculation of [21[ , with a minor correction, to an arbitrary 
number of additional (possibly non-abelian) unbroken charges. 

The observational importance of such an extension lies in the transport terms which 
involve gauge field strengths . For the broken charge, these are inserted as fictitious 
external fields, which serve to increase the power of the entropic argument. In reality, if the 
broken charge is gauged, i.e. in the superconducting case, the gauge fields will be dynami- 
cally excluded from the bulk of the system. However, if we have unbroken gauged charges 
alongside the broken charge, then their associated gauge fields may enter the superfluid, 
rendering these new transport terms observable. In particular, this comment applies to one 
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of the transport terms, which may be called a chiral electric effect: 

J^EE = 0\ c E^U v e p E C G , (1) 

compared to the standard electric conductivity term J^onduct = & ab E£. 

Here c a bc are the transport coefficients, is the normal fluid's four-velocity, ^ is the 
phase gradient of the broken symmetry, which is proportional to the velocity of the superffuid 
part, and is an electric field. We use Greek letters (fj,,u,...) for spacetime indices, and 
Latin letters (a, b, ... ) for charge indices. In (pQ) and later in eqs. fl2}-([B|), we choose to present 
our results with a charge index over the phase gradient This is done in anticipation of an 
extension of the results to the case with multiple broken symmetries, and helps in clarifying 
the index structure of the transport coefficients. We stress, however, that our calculations 
strictly apply only to the case of a single and there is yet much to understand about 
phases with multiple broken charges, even in equilibrium. 

We will present the results for the allowed transport coefficients differently from the 



authors of 



211 ] . First, for concreteness, we will use the transverse fluid frame (for a thorough 



discussion of fluid frames, see {9]). Second, we will show that the results of [2l| can be better 
organized by a different choice of variables, using (s, n a / s, ( 2 ) as a set of independent thermal 
parameters instead of the set (s, fi a /T, ( 2 /T 2 ) used in [21]. Here s is the entropy density, n a 
the normal charge densities, fi a the chemical potentials, T the temperature and ( the length 
of the component of £ M transverse to u 11 [^tJ . Third, we will group the transport coefficients 
in a way which is suggestive of their relation to anomalies. 

Some general remarks are in order here. In conventional hydrodynamics, entropic consid- 
erations only serve to point out which transport terms are allowed. To find the actual form 
of the transport coefficients as functions of the state, one must resort to the microscopic 
theory. In general, transport coefficients are given by n-point functions in the microscopic 
theory, via relations known as Kubo formulas. In this respect, the calculation of was 
seminal: it derived an almost unique form for the chiral transport coefficients in a normal 
fluid, in terms of the coefficient C a f, c of the J J J chiral triangle anomaly. 



In [12j , we noted that the entropic constraints on the normal-fluid chiral terms leave two 



arbitrary constants which were disregarded in llj]. It was later noticed in [2j| that one of 
these constants 7, multiplying a T 3 term in the anomalous current, is ruled out by CPT 
invariance. The other constant (3 a , multiplying a T 2 term in the anomalous current, was 
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evaluated in |24j using Kubo formulas for a theory of free fermions. It was found to be 
proportional to the coefficient of the JTT gravitational triangle anomaly. This result was 
confirmed [25( in a strongly-coupled holographic setup. It may then be suggested that (3 a is 
always related to the gravitational anomaly. Thus, certain transport coefficients may be in 
fact fixed by anomalies, even if the entropic constraints allow them to be more general. 

In the normal fluid, the relationship between the chiral transport coefficients and the 
anomalous 3-point correlators is not entirely clear from the Kubo formula approach. There, 
the coefficients arise from 2-point correlators j^. The 3-point correlator arises effectively 
from the dependence of the fermion propagators on the chemical potentials fi a or the tem- 
perature T. In the entropic approach of [11], the J J J anomaly comes in directly in the 
(non-)conservation law for J£, but the final result for the transport coefficients is removed 
from this by a lengthy calculation. We will propose in section [IV] a new perspective on the 
normal-fluid chiral coefficients, which clarifies their relation to the triangle anomaly and to 
each other. This will be a preliminary step before discussing the superfluid case. 

We will then extend this intuition to the new superfluid transport terms, after grouping 
them in a way which reveals the relevant structure. This leads us to propose a simplified 
form for the chiral constitutive relations at viscous order: 

T$t% = X a A^ v)Xpa u P (a* + a abc &e^ x u p ( ba n XK C + K %^ x u P C ha E cX (2) 

+ B! - £ ( VV/i. + /3"T 2 ) ) (3) 

+ b? c e^u u ( bp ^xC + c abc t^u u Cb P E ca 

^ a „ f^^bcd i nab.. m2 \ , /-a r?u ( 1 r<bcd , , .. , abrrfl 



chiral 



Qc^/Wd + 2/3%T 2 ) + I^fl* Qc 6c Vc// d + /3 fe T 2 ) (4) 



= -y + ^ Qf^^c + 2/3> a r) + B* [^C abc ^ c + . (5) 



The transport terms are presented with the natural generalization to multiple broken 
charges, in order to highlight the index structure of the coefficients. However, we stress 
that terms unique to the case of multiple broken charges are not included. 

In ((21)-©, is the component of £° transverse to w M ; h = e + p is the enthalpy density, 
where p is the pressure and e is the energy density; ir pu is the shear tensor of u p ; w M = 
\e iivpu u v d p u a is the axial vorticity, B% is the magnetic field, and Eg is the combination 



Eg = Eg — TP flu 'V u (fi a /T), where = 5^ + u^u u is the projector orthogonal to u M . z/^ 1 ^ 
is the correction — fi a to the Josephson equation. The transport coefficients x°, a a6c , 
and are arbitrary functions of state, with b^ c + fog 60 satisfying an inequality with 
transport coefficients from the non-chiral sector. If the dynamics is time-reversal invariant, 
we have instead simply bf 10 + b^ = due to the Onsager principle (see section IIII A[) . The 
constant C a b c is the coefficient of the chiral J J J anomaly, and the constant /3 Q is (probably) 
the coefficient of the gravitational JTT anomaly. 

As for the chiral electric conductivity c a b c , we will argue that it's related to the J J J 
anomaly coefficient C a b c . In particular, we suggest that c a b c in (JT]) and (j3J) takes the form: 

abc n dbe ( ra n ° ' \ f cc \ , R s 



A few comments and comparison to 

in 



211 1: The chiral electric conductivity c a bc is denoted 

i n 

2l| as K\x- The coefficients r], k 12 , «2i and K22 from [21] (our x a , bf 30 , b^ c and a abc ) seem to 



be unrelated to anomalies. The coefficient as (will be denoted a a b in this paper) appears to 
be related to a JJT-type anomaly, which does not exist. We therefore expect this coefficient 
to vanish. The coefficient 010 from [2l| (more precisely, 010 — 2(/x/T)o"8 — (C ' /2){ji/T) 2 ) is 
a generalization of the JTT-type anomaly coefficient /3q, where the subscript denotes the 
broken charge. While for a normal fluid the entropic constraints set (3q to a constant, for 
a superfluid they allow it to be an arbitrary function of state. The interpretation in terms 
of the JTT anomaly suggests that the new freedom in the entropic constraints is spurious, 
and /3q is in fact a constant. This conclusion and others are backed by several structural 
arguments, which are presented in section IVl 

The paper is organized as follows. Section HT1 defines our notations and the general frame- 
work of the calculation. In section IHIt we list the transport terms allowed by our calculation 
of the entropic constraints. The calculation is detailed in the Appendix. In section IIVI we 
present our interpretation of the known chiral transport terms in a normal fluid. In section 
IVl we extrapolate from this our educated guesses regarding the interpretation and values of 
the new superfluid transport terms. Section [VI] is devoted to a discussion and outlook. 
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II. FRAMEWORK 



In this section, we outline the framework for the calculation of the chiral transport terms 
allowed by the entropic constraint for a superfluid with a single broken charge and arbitrary 
unbroken charges. The details of the calculation are given in the Appendix, while the results 
are presented in section III. 



A. Definitions and equations of motion 

We consider a superfluid with arbitrary unbroken currents Jf and a single spontaneously 
broken U(l) current Jq. The full set of currents is collectively denoted as The structure 
constants of the charge algebra are f a bc (out of which only may be nonzero). The 
thermal state at each point is determined by the normal velocity the temperature T, the 
phase gradient £ M of the broken symmetry and the chemical potentials /ij of the unbroken 



symmetries. We denote the timelike norm of £ M by £ = \J — In equilibrium, £ M is 
related to the chemical potential fi by the Josephson condition = fi + 0(e), where 
e is a formal small parameter whose powers indicate the number of gradients involved. We 
denote the transverse part of £ M as : 

c„ = p;iu = + w c 2 = w = ti-e , (7) 

where = 6^ + 11^11" is the projector orthogonal to u^. The thermodynamic identities read: 

dp = sdT + n a dfi a + ^-Qd£ 2 (8) 

h = e + p = Ts + fx a n a , (9) 

where p is the pressure, s is the entropy density, fj, a are the chemical potentials, n a are the 
normal charge densities, e is the energy density, h = e + p is the enthalpy density, and Q 
is the (unnormalized) superfluid charge density. The ideal stress tensor, charge current and 
entropy current read: 

T (0)^ = mV + p pHV + Q£»g> ( 10 ) 

jf M = mu" (11) 
4 0) " = n u» - QC (12) 
s (°)" = su» . (13) 
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We take the metric g pv to be curved on the scale of the hydro dynamic gradients, with a 
Riemann tensor R pup a = 0(e 2 ). We couple an external gauge field A p to every current that 
is not already coupled to one. See [3] for the subtleties involved in this procedure in the 
presence of anomalies; as explained there, we use the covariant version of the currents and 
a symmetric anomaly coefficient C a b c . For the broken current Jq , the introduction of the 
external field upgrades £ M = —d^ + from a phase gradient to an arbitrary covector. Its 
curl equals the corresponding field strength: 

F° u = 2d { ^ u] . (14) 

We decompose the field strengths F^ v into electric and magnetic pieces as: 

F« u = 2u { ,E a v] +B; v] El = F^y- B; v = P?P° u F a pa . (15) 

We also define the axial magnetic field vector: 

B a » = l - e ^u v F a pa = l -e^u v B a pa . (16) 

We define a covariant derivative V M which takes into account both the curved metric and 
the gauge fields. We decompose the velocity gradients V pU v into an acceleration a M , a shear 
tensor n pu , a vorticity tensor and an expansion rate V^ii^: 

V^ttx/ = -Uy,a v + 7T /li , + c«v + -V^i/P^ (17) 
a" = u u V u u» (18) 

= p;p;v (pM(T) - l -vyp pv (19) 



00 



P^V [p u a] . (20) 



We also define the axial vorticity vector: 



^ = )f vpa u v d p u a = X -e^u v u pry . (21) 



The ideal equations of motion read: 



v u t^ = f; u j^ + o{e 2 ) (22) 

V M J<°>" = 0{e 2 ) (23) 
u^ = fi + O(e) . (24) 
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The ideal conservation laws ( 122|) -( 123|) can be written as: 



V fl (n u») = V fl (Q^) + O(6 2 ) (25) 

V>^) = 0(e 2 ) (26) 

V„(<m") = 0(e 2 ) (27) 

a" = i - CV^no*")) - ^P^T . (28) 



To obtain eq. ( 12811 for the acceleration, we used the identities (Ej-jS])- This is a slightly 
nonstandard expression, which has some advantages and simplifies our calculation in the 
Appendix. E£ is a combination of the electric field and the chemical potential gradient: 

E% = E% - TP^V^ . (29) 

This is the expression that arises in the standard normal fluid electric conductivity term 
^Conduct = c^Efr. As we will see in section |HB} it plays a role in the entropy constraint. In 
section |V] we will give it an interpretation in the thermal-QFT picture. 
The viscous-order equations of motion read: 

VM 0) " + Ji 1] ") = C ahc E\B^ (31) 
= /i + , (32) 

where is the first-order correction to the stress tensor, is the first-order correction 
to the current, and i/W is the first-order correction to the Josephson equation. C a bc = C(abc) 
is a constant tensor of anomaly coefficients. We work in the transverse frame, defined by: 

Vri^ = 0; V«" = 0. (33) 
B. Entropic constraints on the chiral terms 

Introducing the correction s*- 1 ^ to the entropy current, the second law of thermodynamics 
reads: 

V M (s (0) " + s (1) ^) > . (34) 



S 



We can use eqs. ©-© and (J32J to express V^s^ ^ in terms of V '^ 0> ', V M Ji° )A1 and i/W. 
We can then eliminate V U T^ U and V M< /cs using eqs. (l30]) - (l3Tj) . Eq. is then written 
as a sum of manifestly second-order terms: 



v M ( S (D" + ^ JW" ) > 



T 



(35) 



or, equivalently: 



1. fr^n^ + L T ^u v d vS + j^e; + sz^Vc^ - c abc ^ElB cp 



(36) 



We now wish to find the allowed terms in T^ u , Ja 1 ^, z^ 1 ^ and s 1 - 1 ^ containing e M ^ pcr . For 
normal fluids, the contribution of such chiral terms to the entropy production rate (1351) does 
not mix with the contribution from non-chiral terms. Furthermore, because its sign cannot 



be constrained, the contribution of the chiral terms 
so simple in the superfluid case, as was noted in 



to ( 1351) must vanish. The situation is not 



21] . This is because chiral contributions to 



(|35|) of the form e ,J ' v ' p(7 U fl u l ,CpV (T , where and are some first-order vectors, can mix with 
non-chiral contributions of the form U^U 11 , U^V 1 and V^V P . More specifically, the non-chiral 
contributions may be positive semi-definite with a magnitude that is always greater or equal 
to the magnitude of e^ upa U^u v C, p V a . The coefficient of the chiral contribution can then be 
nonvanishing, without violating the Second Law. The relevant vectors for the role of 
or Vfj, are E°[ (through the J^^E® term in (|36|) ) and 'K ilv C, v (through the T^^ir^ term in 

mi 

III. RESULTS FROM THE ENTROPIC CONSTRAINTS 

We derived the chiral transport terms allowed by the entropic constraint ( 1361) for a su- 
perfluid with a single broken charge and arbitrary unbroken charges. The result reads: 



a A 
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2n a (\ 

- — ( -C bcd Wc fi d + 2a WoT + 2f3 b ^ b T 2 + 7T 3 



(38) 



4r 2 /i C 2 //ig^ a 0flA 2T 2 U a da a d(3 \\ 

h \T d( 2 d( 2 J s \Td(n /s) d(n /s)JJ ' 



(40) 



^2 ^>v r po ■ 



Here x, a, &2, c ab, &a and /3 are arbitrary dimensionless functions of state, while and 
7 are arbitrary dimensionless constants. The partial derivatives with respect to s, n^/s and 
( 2 are taken with (s,n a /s,( 2 ) as the independent thermal parameters; in other words, the 
derivative with respect to n /s is taken at constant s, rii/s and ( 2 , and so on. Note that 
a term of the form e fJlUpeT d u (xu p ( (T ) can be added to (1401) . without changing the dynamics or 
the entropy production rate. 



(140 


Hs 


in 


21| 



nonzero entropy production rate, and will have to be balanced by the non-chiral sector, are 



b a + &| and C[ a b]. The former was noted in 2JJ, while the latter is specific to the case with 



multiple charges. However, we will see shortly that in a time-reversal-invariant theory, the 
Onsager principle imposes the relations: 

b% = c ab = c ba (41) 

This is precisely the condition for which the 6", b\ and c ab terms don't contribute to the 
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entropy production rate. Thus, no counterbalancing from the non-chiral sector is required. 
Let us now demonstrate this relation. 



A. Onsager relations 
In 



21] , it was argued that given time-reversal symmetry, the transport coefficients bf and 
b 2 (their — k u and k 21 ) should be related due to the Onsager principle [26| as k 21 = — k u , 
or, in our terms, = bf. There appears to be a sign error in this relation, as we now show. 

We choose coordinates and a gauge so that locally g^ u = 77^, = (1,0,0,0), £ M = 
(0,^,0,0), and the Christoffel and gauge connection coefficients all vanish. Consider the 
charge q a and the momentum p x along the condensate's velocity The conjugate quantities 
to q a and p x are the electric potential <p a and the velocity v x , respectively. The gradient of 
ip a along the w-axis is d y ip a = —E ay , while the gradient of v x along the z-axis is d z v x = 2tt zx , 
the latter equality holding when d z v x is the only nonzero velocity gradient. The current of 
q a along the y axis is J", while the current of p x along the z axis is T xz . Leaving only the 6" 
and b 2 terms in (l3Tj) - (l38l) (and setting the gradient d y (fi a /T) to zero), we have: 

T xz = ^ble^El = -^b\El = ^bld y <j> a (42) 

r y = c 2 b a 2 ey txz 7T zx = ( 2 b a 2 n zx = ^b a 2 d zVx (43) 

Thus, (( 2 /2)b1 and (( 2 /2)b 2 are mirror-symmetric elements of the kinetic coefficient matrix. 
Now we must pay attention to the time-reversal properties of the relevant quantities. First, 
a time reversal flips the sign of ( x , which is a property of the thermal state. However, in our 
context only the square of ( x enters, so this has no effect. Second, a time reversal flips the 
sign of p x , but not of q a . Therefore, the correct Onsager relation is an antisymmetry of the 
kinetic coefficients, i.e. b 2 = —6". A similar argument applied to the charge currents along 
the y and z axes shows that c a b = Cb a - There, the sign-flip of ( x under time reversal cancels 
with the antisymmetry of e^ upa with respect to the y and z axes. 

IV. REVISITING THE ANOMALOUS NORMAL FLUID 

In this section, we present a new perspective on the known chiral transport terms for a 



normal fluid 
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12[. For the arguments here and in section [V] we regress to the abelian case. 

11 



A careful non-abelian generalization is likely possible, as was done in 12J for the arguments 



of 



Recall the transport terms in the charge current for a normal fluid 121 ] : 

= aX + W ( C abc fi bf i c + 2(3 a T 2 - — i-C bcd ii b ^ d + 2f3 b fi b T 2 + 7T 3 



(44) 



Here C a fc c is the coefficient of the J J J anomaly, (3 a is conjectured to be the coefficient of the 
JTT anomaly, and 7 vanishes due to CPT invar iance. In light of the progression of terms 
C afec -/3 a -7, their charge index structure and the associated factors of fi a and T, we can also 
associate 7 with the pure-gravitational TTT anomaly This interpretation again forces 7 to 
vanish, because such an anomaly doesn't exist in four-dimensional spacetime. 

The combination E% = E% — TP /iI/ V iy (/x a /T) in the electric conductivity term in ( l44"j) 
arises naturally in the context of the Second Law of thermodynamics. This comes about 
through the E^Ja 1 ^ term in the expression ( 135)) for the entropy production rate. The origin 
of the particular combination of vortical and magnetic terms in ( I44p is far less transparent. 
To improve this situation, we propose a certain heuristic way of looking at the transport 
terms. We note that the transport coefficients can be found, via Kubo formulas, from 
correlators in thermal QFT. These can be translated into Euclidean vacuum correlators, with 
a Euclidean metric g^ u , a 1/T periodicity in the imaginary time direction, and an external 
gauge potential [f'u^. However, unlike in the Kubo approach, we will keep discussing non- 
equilibrium quantities such as and B£ directly, instead of translating them into variations 
with respect to external fields. 

With this approach in mind, we expect the system to respond not to the physical gauge 
potential A a ^ but to the hybrid potential A^ = A a ^ + fi a u^. Consider the field strength 
F^ v = 2d[^A^ derived from this potential. Its electric and magnetic parts read: 

e; = P;y = ei - p;d v ^ - (45) 

B% = l - e ^u u F a! „ = B% + 2fi a u» . (46) 

Using the ideal equation ( 1281) in the normal-fluid limit, we find that the electric part ( 1451) is 
in fact proportional to E°: 

E;=U-^)El + 0(e 2 ). (47) 
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Thus, the appearance of in the electric conductivity term is consistent with the approach 
that E® is in fact the field to which the current reacts. 

We should now look for the significance of B£ = B% + 2fi a u fJ- in the chiral terms of 



As a first step towards uncovering it, recall from the Kubo- formula analysis in 23] that the 
chiral terms in J£ are the sum of a "free" piece and a piece multiplied by n a /h: 

J { SZ = J' a " - • (48) 

The point of this decomposition is that the thermal-QFT correlators are directly related to 
J'£ and T" M , rather than to the full combination J^*^. Furthermore, as their names are 
meant to suggest, J'f - and T' M are related to the expectation values (J£) and (T 0fJ- ) of the 
current and the stress tensor, respectively. The T' M term is associated with a reference-frame 
correction between the equilibrium and perturbed states of the fluid. Comparing (148 j) with 
(liljl . we can decompose J'£ and T ,fl into pieces proportional to C a b c , (3 a and 7: 



J a ~ ^abc-J^c) ' > J °- J {(3) 

rp/[_i rplabcpi , n rpla/i , 1 rp/H 

J- — ^abcJ- (C) ~+~ Pa- 1 (/3) +2 M 



(49) 



where the individual coefficients are given by: 

J$ = Sym {p\B°» + ■ Tffi = Sym j/iV (±B^ + l^u^j } ; (50) 

Jg } = 2TV; T$ = T\B a ^ + 4/iV); (51) 

T'fo = 4T 3 . (52) 

Here, "Sym" denotes symmetrization over all charge indices. Consider now the factors ( 150]) 
associated with the J. J J anomaly coefficient C a b c . We can write these as: 

J'$ = Sym f d/i b (B^ + 2/i c u^) = f d^ b B c > = \^ vpa f d(/j (b u u )F^ (53) 
Jo Jo 2 J 

T',^ = Sym / fx a dfx b (B cp + 2//^) = [ ^ a dfi b B c ^ = -e" vpa [ ^ a d{^u v )Ff a . (54) 
Jo Jo % Jo 

The u v inside the integrals is understood to be constant; it is included with //& in the paren- 
theses merely to emphasize the structure of the expression. We understand the contraction 
of (1531) with C a i, c as a relation of the form 5J p /5A b u ~ (7 a 6 c e^ pCT F° commonly encountered 
in the context of anomalies. In this interpretation, /i a w At acts as a constituent of the gauge 
potential A®, as it should when the thermal state is translated into a Euclidean vacuum. 
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The extra factor of \i a in (1341) as compared to (1551) can also be understood to some extent. 
Introducing an external gauge potential affects the stress-energy operator by adding a "po- 
tential energy" term. In particular, the effective gauge potential /i a w M results in an addition 
5T w = ^ S ^j^ NoW; if for 

some reason it is only this addition that enters the chiral 
transport terms, then we have a heuristic explanation for the fact that (I54p is the same as 
(|53|) . with C ahc d^bB^ replaced by C abc fi a d[ibBj?. The integration in eqs. (!53|) - (l54|) can be 
interpreted as the gradual build-up of the relevant thermal-QFT correlators from the state 
with fi a = to the state with fi a ^ 0. This point of view will be utilized in section IVBl 



Recently, the generalization of t 



ie C n hr transport terms was obtained in arbitrary even 
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14 ] have found the generalization of J'^y with [jjj 



spacetime dimensions: two papers 

also giving the generalization of TJ^ CM (for the case of a single U(l) charge). On inspection, 



the results of 



13 



14| satisfy suitably generalized versions of eqs. (j53|) - (l54|) . This lends 



credibility to our emphasis on these relations. For another recent clue regarding the role of 
in the anomalous transport terms, see fl5 |. 

The role of gravitational anomalies in the transport terms is less understood than that 
of the J J J anomalies. Here, we will make do with two modest observations. First, in the 
Euclidean picture, the temperature T has a role with respect to metric variations similar to 
the role of \i a with respect to gauge-potential variations. Indeed, we can consider Euclidean 
spacetimes with a fixed period in the time coordinate, but with different values of the metric 
component goo- Since the inverse temperature 1/T corresponds to the metric length of the 
time period, we then have goo ~ 1/T 2 . This leads to the relation dlngoo = —2d\nT, 
analogous to gL4q = dfi a . Actually, if we wish to consider variations both in the metric 
and in the gauge potential, we should use fi a /T rather than /i a as the quantity associated 
with variations of A a ^. This is because if we vary the metric g^ u in a fixed coordinate system 
without varying A^, what remains constant is not fi a , which is metric-normalized, but rather 
yU a /T, which is the chemical potential's contribution to the Aharonov-Bohm phase § A^dx^ 
over a period of the time coordinate. In section IA 31 of the Appendix, we indeed see that 
(T, Ha/T) is the cleanest choice of variables for deriving the anomalous terms in ((Hj) and 
their superfluid generalization. 

Our second observation is that the /3 a and 7 terms in (jUJ) can be expressed in a form 
tightly related to our expressions ( 153]) -( 154"|) for the C a b c terms. First, let us rewrite ( I53|) -( l54|) 
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in accord with the previous observation regarding the role of /x a /T as opposed to \i a : 

Sym jT5^ j = dJ$; Sym j^TB^} = dlj^ , (55) 

where fi a /T is varied at constant T. Now we point out the following analogous relations for 
the P a and 7 terms in fl5T|) - (l52|) : 



Sym|T^d^} = d(^TS^); (56) 



08) T 

This heuristic observation implies that the transport terms for C a b c , /3 a and 7/2 form a 



Tftd^ = dT'Z . (57) 



sequence, with a missing element between C a ft c and /3 a . This missing element should be 
characterized by a constant rank-2 charge tensor a ao , and its transport terms should read: 

Jg = T5 a ^; 2^*" = H ia TB b) ^ . (58) 

We interpret this sequence as follows. C a bc ~ Tr{C( a C;,G c )} is the coefficient of the J J J 
anomaly, as we know explicitly from the entropic calculation; a ao ~ Tr {G( a Gfy} is the would- 
be coefficient of the nonexisting JJT anomaly (and is indeed absent from (jUjl ); f3 a ~ Tr{G a } 
is the coefficient of the JTT anomaly, as suggested in [24| from a Kubo formula calculation; 
finally, 7/2 is the would-be coefficient of the nonexisting TTT anomaly (which indeed must 
vanish in f )44l) due to CPT invariance). In the above, G a are the generators of the charge 
group in the fermions' representation. 

A generalization of the /3 a and 7 contributions to J'Jf in arbitrary dimensions was re- 
cently discussed in [k| , under the name of "finite-temperature corrections" to the anomalous 
transport terms. A sequence of such terms is found, its length dictated by the spacetime 
dimension. On inspection, this sequence is seen to obey a generalization of eqs. fl56|) - fl57|) . 
Furthermore, the C-term in [l3| is related to the (/3, 7, . . . ) sequence via a generalization of 
( |55|) . again with a "missing link" between the C-term and the rest. We can again interpret 
the sequence in terms of polygon anomalies with a varying number of graviton vertices. The 
slot immediately after the C-term is always empty, because there is never an anomaly with 
a single graviton vertex. 
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V. INTERPRETATION AND EDUCATED GUESSES FOR THE SUPERFLUID 
TRANSPORT TERMS 

In this section, we will use the insights from section |IV]to conjecture a more specific form 
for the transport terms f[5T|) - (l4T)j) . This will bring us to the expressions ©-(jH]) which were 
presented in the Introduction. 



A. Interpreting the a and j3 terms 

Let us return to the superfluid result f )38|) for the chiral part of the current J£. We 
concentrate on the first two terms, involving the vorticity w M and the magnetic field B£. 
We find that these terms reproduce the normal-fluid result (|44p . with two differences. The 
first difference is that the coefficient {3q is no longer a constant, but an arbitrary function of 
state. In the notation of [2l|, it corresponds to a±o — 2(/i/T)cr g — (C/2)(/i/T) 2 . The second 
difference is the introduction of yet another arbitrary function of state, a a . In the notation 
of [2l|, it corresponds to cr 8 . 

To understand better the role of a a , it will be useful to generalize to the case of multiple 
broken charges. We expect that the transport terms for a single broken charge will carry 
through to this more general case, with the trivial addition of a charge index on £ M . In 
addition, new transport terms are likely to appear, involving e.g. factors of & vpa UyC^C^. In 
the present work, we disregard such new terms, and consider only the generalized versions 
of the terms fl37j) - fj40l) . In particular, we find that a a gains a second charge index, becoming 
a a b, for instance, the last term in (T4T)|) becomes (a a b/2)e^ upu QFp a . Following the derivation 
of the vortical and magnetic terms in the Appendix, we see that factors of come to 
replace factors of d^a^. The a a b-contribution to the vortical and magnetic terms in the 
charge current reads: 

jo* = 2 f a W T fy _ ^o)^\ . (59) 

Up to a factor of 2 which can be swallowed into the definition of a^ a b), this reproduces the 
expectation ( 1581) for the missing element in the sequence of coefficients between C a bc and f3 a . 

It may seem at first that since a a and /3o in ( |38l) are arbitrary functions, the whole 
hierarchy of terms from section HVl becomes meaningless. Indeed, why single out the arbitrary 
functions a a and /3 ? We might as well talk about all of (l/3)C abc [i a [ibHc + 2a a /z a /zo^ + 
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2(3 a fi a T 2 + 7T 3 as an arbitrary function. However, a a and /3q really are singled out by the 
entropic calculation. The evidence for this lies in other terms in the constitutive relations 
(I37p - (l39p . which contain partial derivatives of a a and /3q, rather than some other function, 
with respect to the thermal parameters (s, n Q /s, ( 2 ). Furthermore, the derivatives of a a and 
/?o always come together, in the combination B£da a + 2Tu fJ- d(3o. In the spirit of section HV] 
we note that the coefficients of this combination are related by 2Tu fl d(fi a /T) = dB%. This 
reinforces the conclusion that a a and 0o ( OT , more generally, and (3 a ) indeed belong to 
the hierarchy of coefficients described in section IIVt along with C a b c and 7. 

Now, from the entropic calculation we know that C a b c and 7 remain constants in the 
superffuid case. Then the clean hierarchy of coefficients ( 155P - (I58P suggests that at(ab) an d Pa 
are constants as well. By this conjecture, all the partial derivatives of a a and (3 in (]3"7 ]) -( l3"9"]) 
vanish. For /3q, this implies a return to the normal-fluid situation: a constant Po is on a P ar 
with the other constant components of (3 a , which are believed to be the coefficients of the 
JTT anomaly. As for a fab), we can now interpret it as the would-be coefficient of the JJT 
anomaly. Such an anomaly doesn't exist, which leads us to conjecture that in fact ot( a b) — 0. 

For the case of a single broken charge, we are done: otub) collapses back to a a , which 
vanishes by the above argument. For multiple broken charges, we must also consider the 
antisymmetric piece Having come this far, it seems natural to guess that a\ab] vanishes 
as well, so there are no partial-derivative terms in the constitutive relations at all. We note 
in this context that a constant a[ a b\ doesn't affect the dynamics; it only lends an identically 
non-dissipative term to the entropy current. 

We end this subsection by summarizing the constitutive relations with our conjectures 
taken into account: the j3 a are all constants, and a a b vanishes. We invoke CPT invariance 
to remove the 7 term, which has served its rhetorical role. Finally, the transport terms are 
presented with the naive generalization to multiple broken charges, in order to highlight the 
coefficients' index structure; we stress again that terms unique to the case of multiple broken 
charges are not included. The result is given in (j2J)-(|HJ)- 

B. Interpreting the c term: anomalous chiral electric effect 

Let us return to the representation (|53p - (|54p of the anomalous normal- fluid transport 
terms associated with C a b c - Let us consider these expressions in the context of a superffuid 
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with a single broken charge. From the superfruid's point of view, the integration in (l53l) - (|54|) 
builds up from zero the longitudinal part —fiou^ of The result gives us the anoma- 
lous currents in the presence of this longitudinal part. We propose that to obtain the full 
anomalous transport terms for the superffuid, one should proceed analogously to build up 
the transverse component ( p of at constant fi a . This will add the following contributions 
to J' af * and T' M : 

J' a Ne W = ~P^C a0c e^ £ d( u F cpa = -^C a0c e x ^C v F cpa = C a0c e^u u ( p E ca (60) 

1 ft ~ 1 

1 New— o A° 6 / H'a^u^cpa — r,^X U F« e Sv* cpa 

2 Jo 2 (61) 
= C a0c fi a e^u u (pE C(T . 

The projectors P£ are due to the transversality condition (|33|) . The resulting contribution 
to the current reads: 

T (l)ap _ pap, n ltx _ ob I U jl c \ uvpa j- £fi 

J New — J New i 1 New ~ ° I °c h ) a v'spl-'ba ■ 

The relation (jEJ) between E a ^ = E*- P^V v ^ a - /i% and E^ = E*- TP»V u {n a /T) doesn't 
quite hold in the superfluid case, because there is an extra (^-proportional term in the 
acceleration a p . However, a weakened version of the relation does hold: 

P;K = P; (st - ^) Et + 0(e 2 ) , (63) 

where P p = P p — (p( u /( 2 is the projector onto the subspace orthogonal to both and 
We use this to rewrite eq. (162|) as: 

417 = C C ° d (<£ - ! X £ ) (f* ~ ^) e^CA • (64) 
This has precisely the form of the c a b-term from fl38l) . with: 

CoS = C- (s- - - ^tij . (65) 

We propose this as an educated guess for the chiral electric conductivity c a \,. It implies that 
the chiral electric transport term is directly related to the J J J anomaly. Note that our 
expression ( 1651) satisfies the symmetry c ao = Cb a , as required by the Onsager relation ( I4T1) . 

The integration recipe in ( |53|) - (l5l|) and (!60|) - (l6T|) doesn't seem to give the correct results 
if we first build up C M , and then /i a . This makes physical sense: states with fi a = and 
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nonzero are forbidden, since £ M must always be timelike in order to describe the superfluid 
velocity. 

Generalizing our guess f)64p to the case with non-abelian charges and multiple broken gen- 
erators, we get an expression for the superfluid-specific addition to the anomalous transport 
terms: 

4 ] :S = C db£ - (5t - ^ e^u u ( bp E c « , (66) 

where the index b runs over the broken generators. This is equivalent to expression for 
the generalized chiral conductivity c a b c - 



C. Other terms 

The transport coefficients x a i a<lbC ) an d appear to be unrelated to anomalies. We 
expect anomalous terms to be associated with curvatures of the gauge and metric fields. This 
should make them proportional to antisymmetrized derivatives such as F" and, for gravi- 
tational anomalies, uo^ u . On the other hand, the terms corresponding to (x a , a abc , 6" bc , 
are associated with the symmetrized derivative 7T Mi ,. This is not immediately clear for fe* bc , 



but see 



21] or our derivation in the Appendix for its close relationship with b 



abc 
2 ■ 



VI. DISCUSSION 



We analyzed the chiral 
tended the calculation of 



transport terms in relativistic superfluid hydrodynamics and ex- 



2l| to an arbitrary number of additional (possibly non-abelian) 



unbroken charges. We proposed an interpretation of some of the new transport coefficients in 
terms of chiral and gravitational anomalies. We showed that with unbroken gauged charges 
in the system, one can observe a chiral electric conductivity - a current in a perpendicular 
direction to the applied electric field. We proposed an explicit dependence of this conduc- 
tivity on the anomaly coefficient C a b c - Finally, we presented a natural generalization of the 
chiral transport terms to the case of an arbitrary number of spontaneously broken symmetry 
currents. 

There are several open issues for future work. Clearly, our conjectured simplification ()2])- 
(E]) of the transport terms and our proposal for the chiral electric conductivity should be 
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tested with a microscopic calculation. Such a calculation with just one (broken) charge will 
already be a useful check. Also, it will be interesting to have explicit calculations, either 
thermodynamical or microscopic, for the transport terms in a superfluid with several broken 
charges. We expect this more general case to be relevant for nuclear and subnuclear fluids, 
where there are multiple potentially broken generators for the color and flavor symmetries. 

The observational relevance of our results, and indeed of previous results along these 
lines, should be considered. As pointed out in |2l|, the transport terms which aren't related 
to anomalies may have manifestations in nonrelativistic condensed-matter systems. Perhaps 
there is such hope for the anomalous terms as well - though the anomaly is a relativistic 
effect, so is magnetism; nonrelativistic velocities do not necessarily preclude the observation 
of such phenomena. On the particle-physics front, the currently known superfluid phases 
include neutron-star matter and the Color-Flavor locked phase of QCD. The CFL phase is 
particularly interesting for our purposes, since it includes an unbroken gauged generator, 
with respect to which the system behaves as an insulator |8]. This may offer an ideal setting 
for the chiral electric conductivity to be expressed: on one hand, there is an unbroken 
gauge field in the presence of other broken symmetries, as required for the effect; on the 
other hand, since the conventional conductivity vanishes, the chiral conductivity will have 
a strong signature. 
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Appendix A: Derivation of the transport terms from entropic constraints 



1. The first-order terms allowed by symmetries and the ideal equations 

Let us list all the algebraically distinct chiral terms that can appear in the constitutive 
relations: 

i& = up^Cp^ + t a 2 p^c P B p a + wccy + nccc P B p a + + 

(Al) 

+ jfe^u^CpV^ + f^u^C 2 + jte^u u C P E ba (A2) 
"2L = ^iCX + "K»BZ (A3) 

<1) P _ _/fa j{l)an , » I a J3 P _j_ f| Aivpc t TP 

b chiral ~ rp J chiral ' 13 ^ ' °2 1J a ^2 ^ r a P a 

+ s^UvZ&T + ste^u^V^ + s 6 e^u^ p d a C 2 ( A4 ) 

The functions t n , j n , v n and s n are candidate transport coefficients. For the argument that 
the set of terms (1All) - (1A4j) is complete and independent, we refer to [211]. The introduction 
of unbroken charges does not change the reasoning, beyond the trivial addition of charge 
indices. We've omitted from s^ 1 ^ a possible divergence-free term of the form e fJ,vpa d 1/ (xu p ( a ), 
which does not affect the entropic constraints. 

The candidate terms for J^ 1 ^ and s^ 1 ^ are written differently for two reasons. First, 
J^ 1 ^ obeys the transversality constraint (|33|) . while s 1 - 1 ^ does not. Second, the form of the 
Ja 1 ^ terms was chosen for algebraic convenience, while the s 1 - 1 ^ terms were chosen so as to 
ease the calculation of their divergence. 
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2. Contributions to the entropy production rate 



We will now write the contributions to the entropy production rate (136|) arising from (1Alj) - 
(IA4I) . To avoid equivalent terms written in two different ways, we will follow the following 
rules: 

1. We express a p in terms of other vectors, using the ideal equation (J2BJ). 

2. We express factors of V ^ and V M £ M in terms of other scalars, using the ideal equa- 
tions (E5D and m 



3. We express the curl 2d[^ u ] as F® v . 

4. We decompose factors of V ^u u and F" using eqs. (|T7|) and ([15 



5. We avoid the symmetrized derivative V( M ^) whenever possible, using the relations: 

1 

— ( 

2 



6. We eliminate factors of ^^Tfuv, using the relation: 



1 



2\ 



7. We make sure that one index of each e^ upa is always contracted with u^, using identities 
of the form: 

- zx^e^x^UpWa - w x u x e' iVpa Xp,y v ZpU a 
We don't carry out the analogous procedure with £ M , to maintain the standard form 
of the vorticity u p and the magnetic field B%. This leads to a few instances of redun- 
dant forms for equivalent terms; these instances, however, are restricted to the Sj-S\q 
contributions, which will be disqualified independently due to second-derivative terms. 

8. We bring all second derivatives to one of the distinct forms i p V P V p,u V) u p V p F® u and 

9. We note that the gradient of any thermal function can be written as a linear com- 
bination of some k + 2 basic gradients, where k is the number of charges. However, 
we will not bother to do so at this point. Two sets of independent gradients will be 
convenient at different stages of the calculation. 
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10. We note that gradients of the form u^V^x are not independent, since we have 
u^V ^{nil s) = from the ideal equations. At this point, we will leave them as they 
are. 

When commuting derivatives, we take into account the Riemann tensor R^pa and the gauge 
field strength F^ v where necessary. In practice, the gauge field strength never comes up in 
this context. The only place where we must commute gauge-covariant derivatives is in the 
S5 term, where we use the fact that fabcfJ- 0- ^ — 0, since S5 is a group-covariant function of 
fj, a plus scalars. The contributions to the entropy production rate read: 



+ (c (4° - ^k) + *> foe + ^fofo ) 

- | W " foe - fo uUd ^) ~ % B « foe - Ys C ^ duS ) 



+ J ^E ap + J -^B£E afM + ^(^CEav + J -^(,B£eE au - p^d.TuXpK, 

jab . ja jab 

_ JjL_ e ^ P <r V ^ Ui/Cp E aa _ p^d„( 2 u u ( P Ea« - j^e^* E b „u u (pE aa 

- 3 fe^( X 7r x ,u u ( P E aa 

- s\ {ElBZ + £ a ^) + (C,B B V - ^ + c" w ^C) V M 4 

+ S4 ^T (2^ + Bft - 2^C^ - y e^u u C P K) + e^d^u^pdvT 

+ sy,^ fou» + B% - 2u^ - e^uXp (jE b a - ±0 ff r) ) + e^V^C,^ 
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„ " 4^ + 2^) + CX (Sj, + + u^) 



i»E ap (^d a T - -j^E^j + C A (tt A(7 + u; Act )J + V^C^ A V A F ap(T ) 
+ s a 9 U^Cp.Uv (,E ap C^E b ff - ^d a T) + (tt/ + u; p x )B aXa + iu A V A F ap(T 



+ E aix {2^ + B£) - ^Cm^X^*) + e^V p s a 9 u u ( p E a 



+ Sl0 ( 2( p u» ( C ( - ^eA + - ^-u v d v s ) + 2^Cn pu 



+ BffC'fe + + e^%«*(*/ + ^p A ) (Ca (y^ - '^t) + V (ct £ A ) + ^° A ) 
- ££^Cm*UW*CV) + e^fysio^Cpfe + w ctA )C A ■ 

3. Cancellation requirements 

The S7, Sg and Sg terms contain distinct second-derivative terms. Therefore, s-j = Sg = 
Sg = 0. The Riemann-curvature contribution from the Sxq term cannot be canceled with 
anything, so we have sio = 0. 

£5 and ig vanish as the only terms containing ir fll/ ( fJ 'U! u and 'K pv Q l B v a respectively Similarly, 
the tf, tg, tg and ti3 terms cannot be canceled with anything, so they vanish, tu = a and 
t\2 = x don't contribute at all to the entropy production rate. They are therefore arbitrary 
functions of state. If we consider only the parity-odd sector, then the contribution from the 
£" term must cancel with the corresponding contribution from the jg term, giving tf = —jg- 
However, it was noticed in |2l|] that these terms can be counterbalanced by contributions 
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from the parity-even sector. We then have two arbitrary functions tf = b\ and jg = &§■ 
Given time- reversal invariance, the Onsager principle actually gives b% = —6°, removing the 
need for counterbalancing from the parity-even sector; see section IIII Al 

We now turn to the contributions to V^s' 4 without any factors of Collecting the 
coefficients of u^T, u»V ^ a /T), u»d^( 2 , u»E«, B^dJI, B£V^ b /T), B^( 2 and 
we get a generalization of the differential equations in 

9si\ 2si 
T^t — Tfr + 2/x S4 = 0; 



9s 



d(ji a /T)J Tt p 
ds\ 



2Ts a 2 + 2/i s!? = 0; 
+ 2fi s 6 = 0; 



T,Va/T 

2n a s 1 „ 7? 



(A5) 



(77 / ~~ T +5 o s 4-/io (077) 



d(ti b /T)J T>c2 ' u 5 " - nu \d^ b /T) y 
+^6-^of^| =0; 

at > /T^ b /T V at > /T,nb/T 

n°4 + jl_ cabc ^_ _ 2S (a s b) = Q ^ 



h T T 
which can be rewritten as 



\d^ a /T)) Te + T 2 T [ n 

d(s,/T 2 )^ + 2^ Q; 



9 ( 2 /T,ti a /T T 2 

JT = 2T (4 - £ ai ) ; (A9) 

d(s a 2 /T) \ si npf da a \ r abc^. 9 ,(a nb) rAin 

M A* + '° T - T too J T)C2 = C -^7 + 25 a \ (All) 



25 



■ab = C ab Cfic + 2T5^a b) - —Ts\ . 



(A13) 



We renamed S3 = a a } since it will turn out that S3 is an unconstrained function of state. 
Eqs. (jXlQ]) - (TM2j) imply that: 

4 + {sdT + + s 6 dC 2 ) = C^d^ + ^da a + 2*<V>d^ . (A14) 



From this we get: 



~T = 2T2 6 toto + 2 %a>— + P 

Si dT + s a 5 d^r + s 6 d( 2 = -T {^da a + d/3 ) , 



(A15) 
(A16) 



where (3q is an arbitrary function of state and Pi are arbitrary constants. Similarly, eqs. 
(IA6j) - flA8j) now become: 



which can be rewritten as: 



(A17) 



d^-2^(^da« + dp 



C abc ^ c + 48^a b >^ + 2p a )d 



J^2 



(ah) V± 



T 



ito 
T 



d ^ = d ^^C a6c /i a /i 6 /i c + ^a a fi a fi + ^/3> a 
|| = ^C abc fi a fi bf i c + ^a a n a fi + ^P a fi a + 7 



(A18) 
(A19) 
(A20) 



where 7 is another arbitrary integration constant. As noticed in 21], CPT invariance requires 
7 = 0. Summing up and using eqs. fl A9j) and (IA13I) for j a and j% b , we have: 





s~<abc 






s~iabc 


i" 
Jl 


= C abc /ib/i 




2n a , 




h ' 


■ab 
J2 


= C abc fi c 4 



C abc fi a fi b fi c + 2a> a /i + 2/3> a T + 7T 2 



-C bcc Wc/id + 2aWoT + 2(3 b fi b T 2 + 7T 3 
o 



(A21) 
(A22) 

(A23) 
(A24) 



We now collect the contributions to V^s^ of the form e^^U^u^pV^, where and 
are combinations of the independent vectors d^p, d^, V^(/i a /^) and E®. For brevity, we 
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will use the wedge notation U A V, and write gradients as exterior derivatives. On charged 
quantities, we will use the gauge- covariant exterior derivative D, with e.g. D 2 \i a = f a bcF b ii c . 
The cancellation requirement on these terms reads: 



~ fj% n a s A f da a \ 

7 E ' AE ' + (I + T + {w ) l 

,'jf n a s\ ( da a \ \~ u b (ft n a s 6 f da a \ . 




- ^dT A - ^c/T A dC 2 + ds 4 A dT + Ds a 5 A + ds 6 A rf( 2 + A Da a = . 

(A25) 

Considering entropy contributions from the parity-odd sector alone, the first term in (1A25|) 
must vanish, giving = 0, with jg ab ^ unconstrained. However, due to possible mixing 
with the parity-even sector, all of j$ b = c ab can be arbitrary; this is a straightforward 
generalization of the situation with the e^ vpa E^u 1/ ( p tt (T x( x contributions in 21|. On the other 
hand, given time-reversal invariance, c ab is in fact symmetric due to the Onsager principle, 
making the mixing with the parity-even sector irrelevant. The other terms in the first two 
lines of (1A25j) must vanish separately, giving: 



n a S4 


+ 


fda a 




h 


V dT 


n a SQ 




fda a 


+ 




h 


Kd( 2 



/ da a \ 



T,C 2 



(A26) 



Using eq. ( 1A16|) . we write the vanishing of the last line of ( 1A25I) as: 



o = dT a (J^Da a + d Po) + d ( s ^ dT + s I d y + s&dC ' 2 ) ~ s 5 d2 y + TD ^T A Daa 

= dT A (^Da a + dfo) -d(T (y£a a + d/3 ) ) - s a 5 D 2 ^ + TD^A Da a (A27) 
= -fi a D 2 a a - s%D 2 ^ = -fi a f abc F b a c - ^s a 5 f abc F b fi c . 

The f a b c terms vanish, because a a and s§ are functions of /i a plus scalars, and therefore 
commute with fi a under the charge group. Thus, flA27j) is an identity, and doesn't impose 
any further constraints. 

We proceed to examine the remaining contributions to V M s M . Cancellation of the 
Q il uj pj Q J Ev factors requires ft = (fion a /h — 5o)^3- Similarly, the ^B^C^E" factors give 
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it b = (^on a /h — Sfyt^. Collecting the factors of the form ^u^v^V v x in V M s M , we get: 



= u" + ^d,C + ^ (/xoC 2 *s + hv, - 2T Sl ) d, 



no 



2s 4( 9 M T - 2s a 5 V^ - 2s 6 d^ 2 

«" + ;|cU 2 + ^ + - 2T Sl ) 0^ + 2T (^V/ + 9^o)) ■ 

(A28) 

Recalling that u^V ^{rii/s) vanishes due to the ideal equations, we can write this condition 

as: 

= ±ds + ^dC + ^ (/ioC 2 ts + hu! - 2T Sl ) rf^ + + 2T (^Da fl + dft) , 

(A29) 

with arbitrary functions Aj. Choosing (s,n a /s,( 2 ) as the independent thermal parameters 
(and discarding the previous basis (T, fi a /T,( 2 )), we get: 

< A3 °> 

'3 ^ + f ) (A3!) 

2T ^ C 2 + 2T 2 (n a da a d(3 



h h s \T d(n /s) d(n /s) 

\ {^f ahc ii al i hl i c + 2a a fi af x T + 2$y a T 2 + 7 T 3 ^) (A32) 
+ 4T 2 /i C 2 + d/3 \ 2T 2 ( ^ da a , <9/3 



/i V T d( 2 d( 2 J s \T d(n /s) d(n /s) 
Similarly, for factors of the form ^B^u u V u x, we have: 

= T~ s ds + zf dC + Th ^ Chl + hvl ~ Ts2) d! V + A?D 7 + Daa ' (A33) 

with arbitrary A". Again using (s,n a /s,( 2 ) as a set of independent thermal parameters, we 
get: 

da a 

t° = -sT— (A34) 



da a 

tt = ~2T W (A35) 



~ b 2 r l a 



h ' h ' ; s d(n /s) 



h\2 C ^ c + 26 ° a ^ bT + PT ) + ^^W~^d{^/ 
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This concludes the analysis of all the terms in the entropy production rate. Putting every- 
thing together, we arrive at the result fl37|) - fH0|) . 
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